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Subsequent to the early Universe quark-hadron transition the universal baryon number is carried 
by nucleons: neutrons and protons. The total number of nucleons is preserved as the Universe 
expands, but as it cools lighter protons are favored over heavier neutrons reducing the average mass 
per baryon. During primordial nucleosynthesis free nucleons are transformed into bound nuclides, 
primarily helium, and the nuclear binding energies are radiated away, further reducing the average 
mass per baryon. In particular, the reduction in the average mass per baryon resulting from Big 
Bang Nucleosynthesis (BBN) modifies the numerical factor relating the baryon (nucleon) mass and 
number densities. Here the average mass per baryon, me, is tracked from the early Universe to the 
present. The result is used to relate the present ratio of baryons to photons (by number) to the 
present baryon mass density at a level of accuracy commensurate with that of recent cosmological 
data, as well as to estimate the energy released during post-BBN stellar nucleosynthesis. 

I. INTRODUCTION 

Big Bang Nucleosynthesis (BBN) is regulated by the competition among various two-body collisions leading to 
the creation of complex nuclides (mainly helium-4, with trace amounts of deuterium, helium-3, and lithium-7) built 
from neutrons and protons. The nucleon number density (ub) is key to the BBN-predicted primordial abundances. 
In contrast, for studying the growth and evolution of structure in the Universe, the baryon (nucleon) mass density 
(pb), including the small contribution from the accompanying electrons, plays a more direct role. Of course, the two 
parameters are directly related through the average mass per baryon: pb/"-b = "t-b- Given the high precision of recent 
cosmological data, it is important to match the adopted value of the mass per baryon to the accuracy demanded by 
the data when relating the baryon-to-photon ratio to the baryon mass density parameter. It is shown here that as 
the Universe expands and cools, me evolves from a value of me = 939.17 MeV at very early times, shortly after 
the completion of the quark-hadron transition, to a value of tob — 938.88 MeV just prior to the onset of BBN at 
T w 70 keV, a few minutes after the expansion began. The post-BBN average mass per baryon has decreased to 
niB = 937.12 MeV, leading to r]io = 10^''(nB/n^)o = (273.9 ± 0.3) Ob^i^, where ^Ib is the present ratio of the baryon 
mass density to the critical density and h is the present value of the Hubble parameter in units of 100 kms~^Mpc^^. 

II. THE AVERAGE MASS PER BARYON 
A. Pre-BBN 

Prior to Big Bang Nucleosynthesis, at temperatures below the quark-hadron transition, the baryons present in the 
Universe are nucleons: neutrons and protons (m„ — 939.565 MeV, nip — 938.272 MeV [1, 2]). In weak equilibrium 
at high temperatures there are nearly as many neutrons as protons: n/p — exp(— Am/T) —^ 1 for T ^ Am = 
1.2933 MeV [1]. In this unachievable limit, jtib {ran + ™h)/2 = 939.17 MeV. Note that the hydrogen mass 
(mn = 938.783 MeV [2]) is used here since in the later evolution of the Universe the proton and electron go hand in 
hand. However, even at T = 100 MeV, the ratios of baryons to protons (= hydrogen = H) by number and by mass 
are 

T = 100 MeV : ne/nH = l + n/p = 1.98715, pb/pb = 1 + 1.000833(n/p) = 1.98797, (1) 

so that the average mass per baryon remains at tub = 939.17 MeV. As the Universe cools further, the ratio of neutrons 
to protons shifts even more in favor of the lighter proton, reducing the average mass per baryon. At temperatures 
below ^ 1 MeV the neutron to proton ratio deviates from (exceeds) its equilibrium values, and at the onset of BBN 
at a temperature of ~ 70 keV, detailed calculations [3] reveal that n/p k 1/7. At this time, 

r = 70kcV: nB/nH = l + n/p« 8/7 = 1.14286, pb/ph = 1 + 1.000833(n/p) « 1.14298, (2) 

so that the average mass per baryon has now decreased to niB ~ 938.88 MeV. 
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B. Post-BBN 

During BBN the free neutrons and protons are rapidly transformed into light nuclides, resulting in a mixture 
dominated by hydrogen and hclium-4, along with trace amounts of deuterium and hclium-3 (any produced during 
BBN decays to ^He). The mass-7 nuclides (^Li, ''Be) are the only others produced in astrophysically interesting 
abundances, but the BBN-predicted abundance of mass-7 is so small as to ensure its contribution to the analysis here 
is negligible (for a review of BBN and further references, see Stcigman (2005) [4]). Since the binding energy of these 
light nuclides has been radiated away, the average mass per baryon has decreased. Until stellar processing begins, 
much later in the evolution of the Universe, the baryon mass and number densities are dominated by H, D, ^He, and 
''He, and the average mass per baryon is frozen at its post-BBN value. 

Introducing the abundance ratios by number with respect to hydrogen, yi = ni/nn, the post-BBN baryon number 
and mass densities, may be related to the corresponding number and mass densities of hydrogen, 

nB/nH = l + 2?/2 + 3j/3 + 42/4 + ... , (3) 

Pb/Pk = 1 + {m2/mu)y2 + {m3/mu)y3 + {mi/mn)yi + ... • (4) 

Since any electrons which remain after annihilation are coupled clectromagnetically to the charged nuclei, the 
appropriate masses to be used in eq. 4 are the atomic masses of the neutral atoms. The neutral atom masses adopted 
here from Audi et al. (2003) [2] are: m2 = 1876.12 MeV, tos = 2809.41 MeV, rrn = 3728.40 MeV. 

As a zeroth approximation, ignore the contributions from D, ^Hc (and ^Li), and assume that the primordial ^He 
abundance is j/4 = 1/12 (so that the primordial **He mass fraction, defined as Yp = 4^4/(1 + 4^4), is Yp = 0.250). 
This assumes that all neutrons available at BBN are incorporated into *He and corresponds to a neutron to proton 
ratio, evaluated just prior to the onset of BBN, n/p = 1/7. In this approximation, ne/nn = 4/3 and Pb/ph = 1.3310, 
so that niB = 937.11 MeV. For arbitrary Yp (see Cyburt (2004) [5] and Serpico et al. (2004) [6]), 

mB/mH = 1 - (1 - i(!I^))Yp = 1 - 0.007119Yp, (5) 

4 JTlH 

or 

TUB = 938.112 - 6.683(Yp - 0.250) MeV. (6) 

For Yp = 0.250, this yields me = 937.11 MeV, in perfect agreement with the result above. 

It is possible to improve upon this approximation by utilizing the BBN-predicted, primordial abundances of the light 
nuclides. For example, the WMAP-Only, 3-year data from Spergel et al. [7] correspond (see §111 below) to a baryon- 
to-photon ratio of ryio = 6.12^0 25. Using this estimate in the fitting formulae for the BBN-predicted abundances 
from Kneller & Steigman [3], the primordial abundances of the light nuclides may be calculated. For 7710 = 6.12, the 
BBN-predicted abundances are: y2 = 2.56 x 10"^, j/s = 1.05 x 10-^ and Yp = 4y4/(l + 4?y4) = 0.2482. With these 
abundances, Ub/tih = 1.3302 (in contrast to 4/3 in the zeroth approximation above) and ps/pa = 1.3279 (compared 
to 1.3310 above). However, their ratio, the average mass per baryon, changes from the zeroth approximation by only 
0.1% to rriB = 937.12 MeV. It should be noted that this small difference is entirely driven by the slightly different ^He 
abundance (0.2482 vs. 0.250) and is completely consistent with eqs. 5 & 6; in the ratio of pB to ub, the contributions 
from D and ^He cancel to a few parts in 10^. Thus, from shortly after the quark-hadron transition (T « 100 MeV) 
until stars form, reinitiating nucleosynthesis, the average mass per baryon has decreased by 2.05 MeV, from 939.17 
MeV to 937.12 MeV. 



C. Post-Dark Ages 

During the very recent evolution of the Universe, initially small fluctuations in the mean matter density have 
been amplified, resulting in nonlinear perturbations which collapse under their own gravity, separating such regions 
from the average expansion of the Universe and leading to the formation of stars, galaxies, and cluster of galaxies. 
Prior to the gas being cycled through stars where hydrogen is burned to helium and beyond, the mean mass per 
baryon is preserved at its precoUapse/post-BBN value. However, as stellar nucleosynthesis proceeds, baryons are 
incorporated into ever more tightly bound nuclei whose binding energy is radiated away. As a result, locally, the 
average mass per baryon decreases, and ttib is now distributed inhomogeneously throughout the Universe. So, while 
the spatially averaged value of me will be smaller than the post-BBN value, local values of tob will depend on 
local histories of stellar nucleosynthesis. As an illustration we may adopt the recently revised solar abundances from 
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Asplund, Grevcsse and Sauval (2005) [8] (AGS) {Xq = 0.7392, Yq = 0.2486, Zq = 0.0122). The ratio of baryons to 
hydrogen by mass is (Mb/Mh)© = Xq^. For these values, 2/40 = 0.08468 (4^/40 = 0.3387). Adopting the Gciss and 
Gloeckler (1998) [9] proto-solar D and '^He abundances, (1 + 2y2 + 3?/3 + 4^/4)0 = 1.3388. Accounting for the baryons 
incorporated in all the heavy elements increases the ratio of baryons to hydrogen by number to (A^b/A^h)© = 1-3554, 
leading to (me/mH)© = 0.9981 and a local value of me© = 936.99 MeV. A word of caution is called for here. 
The AGS abundances [8] are photospheric and certainly not representative of the proto-solar abundances. A better 
choice may be the BP04-I- presolar abundances adopted by Bahcall, Serenelli and Pinsonneault (2004) [10] (BSP) 
(A© = 0.71564, Yq = 0.26960, Zq = 0.01476). In this case, 2/4© = 0.09486 (4y4© = 0.37943); since (Z/A)©(BSP) 
= 1.25(Z/A)©(AGS), the heavy element abundances (by number with respect to hydrogen) are scaled here by this 
factor. For BP04+, (A^b/A^h)© = 1.40021, leading to (mB/wH)© = 0.99796 and rriB© = 936.87 MeV. Thus, even in 
the chemically enriched solar vicinity, the average mass per baryon has been reduced from its value during the dark 
ages by only ^ 0.26 MeV. 

For a BP04+ mix of abundances, along with the assumption that the helium abundance increases in proportion to 
the increase in metallicity {AY oc AZ = Z), it follows that the average mass per baryon decreases with increasing 
metallicity as Atub ~ 0.2Q{Z/Zq) MeV. This decrease in the average mass per baryon mirrors the energy released 
when hydrogen is burned to helium and beyond. With the simpliiying assumption that stellar nucleosynthesis occurs 
at redshift 2;,, and where Z is the large-scale average metallicity, the corresponding present day ratio of the density 
of the energy released to the baryon mass density is 

AOb _ AmB 1 _ 2.7 X 10-'*(Z/Z©) 

~ ruB (l + z*) (1 + z,) ■ ^' 

For the Spergel et al. [7] value of Qb/i^ = 0.02233, AQb/i^ « Q.lxlQ-^{Z / Zq){1 + z^)'^; for h = 0.7, this corresponds 
to Af^B ~ 1-2 X 1Q~^{Z/Zq){1 + z*)~^. It is interesting to compare this to an estimate of the observed energy density 
in background light. According to Fukugita and Peebles (2004) [11], the sum of the optical and far-infrared (FIR) 
backgrounds is ^opt/fir ~ 2.4 x 10~^, so that (see Hauser and Dwck (2001) [12]) 

AOb ^ HZ/Zq) ^ 

^OPT/FIR (1 + Z*) 

suggesting that stellar nucleosynthesis should occur at redshifts 1+ z», ^ ^{Z/Zq) or, that Z/Zq ^ 0.2(1 -h 2*), if the 
observed background is taken as a lower limit to AQb- 



III. THE BARYON NUMBER AND BARYON MASS DENSITY PARAMETERS 

As the Universe expands in the post-BBN era, the number of baryons in a comoving volume is conserved and the 
average mass per baryon is frozen at its post-BBN value until stellar nucleosynthesis begins. Since the baryon number 
and mass densities decrease as the Universe expands, it is convenient to express the baryon number and mass densities 
in terms of parameters which remain unchanged as the Universe expands. The baryon mass density, pB, is usually 
written as a fraction of the critical mass density, defined by pc = SHq/SttGn, Ob = Pb/Pc- Hq is the present value of 
the Hubble expansion rate parameter, expressed as _ffo = 100ft kms~-'^Mpc~-'^. According to the Particle Data Group 
(PDG) [1], Pc = 1.0537 X IQ^^h^ GeVcm^''. The fractional error in pc/h'^^ 1-5 x 10^^, is entirely due to the uncertainty 
in the experimental vahic of Newton's gravitational constant (Gn = 6.6742 ± 0.0010 x 10^^^ m'^kg~^s~^ [1]). This 
uncertainty, while not entirely negligible, is subdominant and will be ignored here. The ratio of the present baryon 
mass density, pB, to the hydrogen mass (toh = 0.938783 GeV), is related to Ob/i^ by 

Pb/toh = 1-1224 X 10"^ VLbH'^ cm'^ (9) 

It is conventional to compare the baryon number density, tib, to the number density of the Cosmic Background 
Radiation (CBR) photons, n^. To high accuracy this ratio, r}\Q = 10^°(nB/n-y)o, is unchanged from the end of 

BBN to the present. Adopting the present CBR temperature from Mather et al. (1999) [13], T^o = 2.725 ± O.OOIK 
= 2.725(1 ± 0.0004)K, the corresponding present epoch number density of CBR photons is 

n^o = 410.50(T^o/2.725K)-'' cm--'' = 410.50(1 ± 0.0011) cm'^ = 410.5O ± 0.45 cm"^. (10) 

Note that for the same adopted temperature the value listed in [1] (410.4), is wrong and has been corrected [14]. 
Combining the results above, rjio may be related to Ob/i^, 

t^io/Ob/i^ = 273.42(TOB/mH)"^(2.725K/T^o)^ = 273.42(1 - 0.007119Yp)-i(2.725K/T^o)^- (H) 
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This result is consistent with eq. 4.19 of Serpico et al. (2004) [6]; with the updated value of Gn [1] their coefficient 
273.49 becomes 273.42. For T^o = 2.725 ± O.OOIK, 

Vw/^Bh"^ = 273.9 ± 0.3 + 1.95(Yp - 0.25). (12) 

This numerical conversion factor (for Yp = 0.25), ryio/OB^i^ = 273.9 ± 0.3, is a key result of this Letter. Note that 
even for the solar value of the average mass per baryon derived above using the BSP [10] abundances (BP04+), this 
factor increases only very slightly, from 273.9 to a local value of 274.0. 

The ratio of 7710 to ilsh^ can be calculated to an accuracy of 0.1%, considerably better than the current 3 - 4% 
accuracy of the flBh^ determination from WMAP. However, to avoid contributing an avoidable error to 7710 as inferred 
from non-BBN data, this coefficient, 273.9, should be employed when converting from UbH^ to 7710. For comparison, 
while Trotta and Hansen (2004) [15] use CAMB/RECFast for their numerical results, in their analytic results they 
employ a conversion coefficient closer to 275, which differs from the result here by 0.4%. The Ichikawa and Takahashi 
(2006) [16] conversion coefficient of 273.49 (based on Serpico et al. [6], using an older value of Gn) is closer to the 
correct result, but still differs from it by more than the uncertainty of 0.3. For the WMAP-Only, 3 year data, Spergel 
et al. (2006) [7] quote (in Table 5) Ob/i^ = 0.02233t!^;o^ogi, corresponding to mo = 6.116t!^;249 « 6.12^°^° (not to 
the value. t/io = 6.0965 ± 0.2055, listed in Table 4 of reference [7] which corresponds instead to ryio/f^B^^ = 273.0). It 
is this value, rjio = 6.12, which was used in §nB to calculate the primordial abundances needed for a more accurate 
computation of the post-BBN average mass per baryon niB. 

Finally, it should be noted that the confrontation of cosmological models with the CMB data depends on both the 
baryon number density (through the number density of electrons) and the baryon mass density. As a result, such 
analyses depend on the numerical value of the coefficient relating the two. The codes most commonly employed to 
analyze the CMB data, CAMB [17] and CMBFast [18], use RECFAST [19]. The latter assumes that toho = 4toh, so 
that mB = m-H and, it adopts a numerical value for mu which is too large by 0.108 MeV. In addition, the codes employ a 
value for Newton's constant, Gn = 6.67259 x 10~^^ m^kg~^s~^, which differs slightly from the currently recommended 
best value [1, 14]. As a result, the conversion factor built into these codes corresponds to 7710/^^3/1^ = 273.46 (provided 
that T^o = 2.725K is adopted), which differs from the more precise result, 273.91, by 1.5 times the error contributed 
by the uncertainty in T^o. It should be noted that since Spergel et al. [7] use CAMB [17] (and RECFAST [19]) in 
their analysis, in principle their derived value of QbIi^ = 0.02233 corresponds to ijio = 6.106 « 6.11. However, this 
small difference in rjio doesn't change the BBN-predicted primordial *He abundance (to four significant figures) and, 
as a result, has no effect on the conversion factor calculated here. 
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